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Abstract 

The double-radiative process e + e~ — > Z77 — > qq77 where the two hard photons 
escape detection at low polar angles into opposite directions, is studied in 0.62 fb _1 
of data collected with the L3 detector at LEP at centre-of-mass energies between 
188.6 GeV and 209.2 GeV. The cross sections are measured and found to be con- 
sistent with the Standard Model expectations. 
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1 Introduction 



One of the most copious sources of events in e + e~ collisions at LEP above the Z resonance 
is the process e + e~ — > qq, with a cross section of about 100 pb. The effective centre-of-mass 
energy, y/s', at which this hadron production takes place does not necessarily correspond to 
the centre-of-mass energy of the LEP machine, y/s, owing to the emission of one or more hard 
initial-state-radiation (ISR) photons by the incoming electrons or positrons. These photons are 
most likely emitted along the beam line, in the low polar-angle regions of the detectors which 
are not instrumented and, therefore, escape detection. 

The cross section of the e + e~ — ► qq process was measured [1,2] and found to be in agreement 
with the Standard Model predictions for both a subsample of events with values of xfs 1 close 
to y/s and a more inclusive sample extending to lower values of V^- The emission of ISR 
photons often implies x/P ~ mz, where mi = 91.19 GeV is the mass of the Z boson. This 
phenomenon is commonly called "radiative return to the Z " . The process e + e~ — > Z7 — > qq7, 
where a hard ISR photon is responsible for such radiative return to the Z, was studied in 
detail [4,5]. Events in which the photon was visible in the detector, were used to constrain 
possible anomalous triple-couplings between neutral gauge bosons [4]. Events with either a 
detected photon or a low-angle undetected photon were used to reconstruct the mass of the Z 
boson and validate the analysis tools used in the measurement of the mass of the W boson [5]. 
The e + e~ — > Z77 — > qq77 process, where both ISR photons were visible in the detector, was 
first observed by the L3 collaboration [6]. The cross section of this process was then measured 
for y/s = 130 — 209 GeV and found to be in agreement with the predictions [7] . 

This Letter extends the study of the e + e~ — » Z77 — > qq77 process to the case in which 
both ISR photons are emitted at low polar angles and are therefore not detected. In particular, 
the case is considered in which the two photons are emitted on opposite sides of the detector, 
with comparable transverse momenta. In this topology, the two jets originating from the Z- 
boson decay are back-to-back. In the following, this process is denoted as "double-radiative 
return to the Z " . This study complements the previous studies of the e + e~ — > qq process in 
a very specific phase-space region and allows further tests of Monte Carlo simulations of ISR 
photons in hadronic events. Moreover, final states with two back-to-back hadronic jets and 
missing energy are a signature of the near-threshold production of the Standard Model Higgs 
boson, H, at LEP in the reaction e + e~ — > ZH. In this case, the missing energy is due to a Z 
boson decaying into neutrinos and the jets to the Higgs boson. In addition, manifestations of 
New Physics in the production of an invisibly-decaying Higgs boson in association with a Z 
boson decaying into hadrons would also give rise to the same final state. Finally, similar event 
topologies are predicted by Supersymmetry. Therefore, a study of double-radiative return to 
the Z with unobserved photons validates the background Monte Carlo simulations for those 
searches. 

This analysis selects e + e~ — > qq events with two or more hard ISR photons satisfying the 
following phase-space criteria: 

E 112 > 5 GeV (1) 
|cos0 7l 2 | > 0.96 (2) 

Iv^-mzl < 2r z (3) 
cos# 7l x cos# 72 < (4) 
< 0.1v/i, (5) 
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where E li , 6 li and p^. are the energy, polar angle and momentum in the plane transverse to the 
beams of the photon i, respectively. Tz denotes the width of the Z boson, 2.49 GeV [3]. If more 
than two ISR photons are present in the event, this signal definition is applied to the two most 
energetic ones. These criteria select about 6% of the phase space of the e + e~ — ► qq process, 
corresponding to a cross section of about 5.5 pb in the y/s range explored at LEP. Figure [T] 
illustrates the complementarity of this phase space with those covered by the analyses of the 
e + e~ — > qq, e + e~ — Z^y — >■ qq7 and e + e~ — > Z77 — > qqj^f processes described in References 1, 4 
and 7. 

2 Data and Monte Carlo Samples 

This measurement is based on 0.62 fb _1 of data collected with the L3 detector [8] at LEP 
in the years from 1998 through 2000 at centre-of-mass energies between y/s = 188.6 GeV and 
y/s = 209.2 GeV, as detailed in Tabled In the last year of data taking, the LEP centre-of-mass 
energy was routinely increased while the beams were colliding in order to enhance the sensitivity 
of the search for the Standard Model Higgs boson, exploring the range y/s = 202.5 — 209.2 GeV. 
In the following, this last data sample is split into two energy ranges. 

The KK2f [9] Monte Carlo program is used, with default options, to generate a total of 1.9 
million e + e~ — > qq events which can contain one or more hard ISR photons, at the centre-of- 
mass energies listed in Table ^ These events correspond to about 35 times the luminosity of 
the data and cover a phase space much larger than that of the criteria (1) — (5). If at least two 
ISR photons which satisfy the criteria (1) — (5) are present in an event this is treated as signal, 
otherwise it is considered as background. This distinction between signal and background is 
performed on generated variables, before any event simulation and any application of detector 
resolutions. 

Other background processes are generated with the Monte Carlo programs PYTHIA [10] for 
e + e~ — ► Ze + e~ and e + e~ — > ZZ, KK2f for e + e _ — > r + r~, PHOJET [11] for hadron production 
in two-photon collisions and KORALW [12] for W-boson pair production except for ez/ e qq' final 
states, generated with EXCALIBUR [13]. The hadronisation process for signal and background 
events is modelled with the JETSET [10] program. 

The L3 detector response is simulated using the GEANT [14] and GHEISHA [15] programs, 
which model the effects of energy loss, multiple scattering and showering in the detector. Time- 
dependent detector efficiencies, as monitored during data-taking periods, are also simulated. 

3 Event Selection 

The event selection proceeds from a sample of high-multiplicity events. Events containing 
photons, electrons or muons with energies above 20 GeV are removed in order to reduce the 
backgrounds from e + e~ — > qq events with ISR photons in the detector and events containing 
W bosons which decay into leptons. The visible mass, M vis , and the visible energy, E vis , of 
these events are required to satisfy 50 GeV < M vis < 140 GeV and 0.4 < E vis / y/s < 0.65, to 
reduce both e + e~ — > qq events without missing energy due to ISR photons and most events 
from two-photon collisions. The latter cut is illustrated in Figure EK- Events from two-photon 
collisions are further suppressed by requiring | cos thrust I < 0.96, where thrust is the angle 
between the thrust axis and the beam line. Events are then reconstructed into two jets by 
means of the DURHAM algorithm [16] and the signal signature of two back-to-back jets is 
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enforced by requiring the angle between the two jets, 6^, to satisfy 6^ > 1.5 rad. Finally, the 
sum of the momenta of the two jets in the plane transverse to the beams, pt, must be less 
than 0.2E v i S . This cut, shown in Figure Eb, accounts for the fact that all missing momentum 
in signal events is due to the two ISR photons nearly collinear with the beam particles and 
therefore directed along the beam line. After this pre-selection, 17208 events are selected in 
data, well consistent with the 17151 events expected from Monte Carlo simulations, of which 
13% are from signal and 87% from background. The background is almost entirely composed 
by e + e~ — > qq events which do not satisfy the signal definition (1) — (5). Small contributions 
arise from four-fermion production and hadron production in two-photon collisions. The signal 
efficiency at this stage of the analysis is 68%. 

Three additional cuts are devised to reduce the residual background and enhance the signal 
component in this sample. The energy of the most energetic jet must be greater than 0.4^; 
the angle between the two jets in the plane transverse to the beams, 9^, is required to be 
0jj > 2.9 rad, as shown in Figure Eb; the polar angle of the jet closest to the beam line, 
should be such that I cos 0j^| < 0.85. Finally, two of the pre-selection criteria are tightened: 
9jj > 1.95 rad and 70 GeV < M vis < 110 GeV, as shown in Figures Eli and El respectively. The 
former criterion is extremely efficient in removing the background from one-photon radiative 
return to the Z boson, which is characterised by a larger boost than the signal and therefore a 
smaller jet opening-angle. 

After these selection criteria, 1672 events are selected in data while 1684 are expected from 
Monte Carlo simulations, of which 61% are from signal, and 39% from background, as detailed 
in Table El Three quarters of the background are due to e + e~ — > qq events which do not satisfy 
the signal definition (1) — (5). The remaining background is due to four-fermion production and 
hadron production in two-photon collisions. The average signal efficiency is 31%. 

The distribution of M vis , shown in Figure El presents a clear enhancement at m<z, as expected 
for signal events. 

4 Systematic Uncertainties 

Several sources of possible systematic uncertainties are considered, and their effects are sum- 
marised in Table 01 

Monte Carlo simulations might not perfectly reproduce the tails of the variables used in the 
event selection owing to, for instance, non-linearity in the modelling of the calorimeter response 
or a bias in the determination of jet directions close to the edge of fiducial volumes. To assess 
this effect, the analysis is repeated by removing one selection criterion at a time. In addition, 
a 0.5% uncertainty in the jet energy-scale and a 2% uncertainty in the determination of the jet 
angles are also considered. 

The signal and background events from the e + e~ — > qq process are generated taking into 
account the interference between ISR photons and those emitted in the final state. The anal- 
ysis is repeated by using a Monte Carlo sample without this interference and the difference 
with the original result is used as an extreme systematic uncertainty on the modelling of this 
phenomenon. 

The cross sections are measured by assuming a fixed background level, as discussed below. 
Uncertainties in the background cross sections are therefore a possible source of systematic 
uncertainty, which is estimated by repeating the analysis with a variation of 10% for the cross 
section of the e + e~ — > ez/ e qq' process, 5% for e + e~ — > qq events classified as background, 5% 
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for the e + e~ — > ZZ process, 5% for the e + e~ — > Ze + e~ process and 0.5% for W-boson pair 
production. 

Finally, statistical uncertainties related to the limited amount of Monte Carlo events used 
to describe the signal and the background processes are included as systematic uncertainties. 
The total systematic uncertainty on the signal cross section varies between 5.3% and 7.7%, 
depending on the centre-of-mass energy. 

5 Results 

The signal cross sections are determined for each centre-of-mass energy by fitting the observed 
distributions of M v i S . Two components are considered, both with a shape fixed to Monte Carlo 
expectations: a signal component with a free normalisation, and a background component with 
fixed normalisation. The results are listed in Tableland plotted in Figure together with the 
corresponding statistical and systematic uncertainties. A good agreement with the predictions 
of the KK2f Monte Carlo, also given in Table El and Figure EJ is observed. These predictions 
have an uncertainty of 3%, which includes a statistical component and the uncertainty from 
higher-order corrections, estimated following the suggestions in Reference 9. 

To further compare the results and the expectations, the ratio between the measured, er, 
and the expected, cr t h, values of the cross section is calculated for each centre-of-mass energy. 
These values are then averaged, by assuming all systematic uncertainties to be fully correlated, 
with the exception of those due to the limited Monte Carlo statistics. The result is: 

a/a th = 0.98 ± 0.04 ± 0.06 , 

where the first uncertainty is statistical and the second systematic. 

In conclusion, the cross section of the process e + e~ — ^ Z^y — >■ qq77, where the two photons 
are emitted in the phase space defined by the criteria (l)-(5), is measured with an accuracy of 
7% and is well reproduced by the current simulations of ISR in hadronic events. This finding 
validates the estimate of the background from events with two back-to-back jets with mass 
close to the mass of the Z boson both in the searches for Higgs bosons of the Standard Model 
and beyond and for other manifestations of New Physics. 
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-y o ^uc v j 


188 6 


191 6 


195 fi 


199 5 


201 7 


202 5 — 205 5 


205 5 — 209 2 


C (pb-i) 


176.0 


29.5 


83.4 


81.4 


36.7 


77.5 


138.6 



Table 1: Centre-of-mass energies and corresponding integrated luminosities, C, con- 
sidered in this analysis. The last two energy ranges correspond to the average 
centre-of-mass energy values < \/s >= 204.8 GeV and < \/s >= 206.6 GeV, re- 
spectively. 





A^Data 


Nmc 


^Sign 


^Back 




a 


0"th 


(GeV) 










(%) 


(pb) 


(pb) 


188.6 


632 


660.1 ±6.9 


372.3 ±4.6 


287.9 ±5.1 


36.1 ±0.4 


5.31 ±0.39 ±0.37 


5.9 


191.6 


84 


97.3 ±2.1 


56.7 ± 1.5 


40.6 ± 1.5 


34.7 ±0.9 


3.63 ±0.88 ±0.44 


5.7 


195.5 


238 


236.5 ±4.2 


149.7 ±2.8 


86.8 ±3.1 


32.3 ±0.6 


5.51 ±0.55 ±0.32 


5.5 


199.5 


216 


196.3 ±3.9 


126.6 ±2.5 


69.6 ±3.0 


29.6 ±0.6 


6.05 ±0.60 ±0.31 


5.3 


201.7 


82 


80.9 ±2.1 


52.8 ± 1.6 


28.1 ± 1.4 


27.8 ±0.8 


5.32 ±0.87 ±0.33 


5.2 


204.8 


147 


153.2 ± 1.6 


102.2 ± 1.0 


50.9 ± 1.2 


26.4 ±0.3 


4.81 ±0.58 ±0.27 


5.1 


206.6 


273 


259.9 ±2.3 


173.5 ± 1.3 


86.5 ±2.0 


25.5 ±0.2 


5.33 ±0.46 ±0.27 


5.0 



Table 2: Number of events observed in data after the final selection, A^ata, compared with the 
total number of events expected from Monte Carlo, N MC . The number of signal events expected 
from the KK2f Monte Carlo, iVsi gn , is also given, together with the number of background 
events, iV Bac k- The selection efficiency, e, is also listed, together with the measured, a, and 
expected, <T t h, signal cross sections. The uncertainties on Nmc, Nsi gn , A^ack and e correspond 
to the statistical uncertainty of the Monte Carlo. The first uncertainty on a is statistical, the 
second systematic. 
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Source 


Effect (%) 


Selection criteria 


2.3 


Jet energy scale 


3.0 


Jet angle 


3.0 


ISR/FSR interference 


1.2 


Background normalisation 


1.4-4.4 


Monte Carlo statistics 


1.2-3.9 


Total 


5.3-7.7 



Table 3: Systematic uncertainties on the signal cross section. 
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Figure 1: Distributions at generator level of the absolute difference of the two photon 
polar angles versus the energy of the most energetic photon for a) the signal, b) the 
e + e~ — > qq process for yfs* / yfs > 0.85 and yfs' > 60 GeV, c) the e + e~ — > qq process 
for y/s' > 60 GeV, d) the e + e~ — > Z7 — > qq7 process, and e) the e + e~ — > Z77 — > 
qq77 process. Only events with at least two photons with energies greater than 
1 MeV are shown, for a sample at y/s = 189 GeV. The histograms show the fraction 
of the cross section of each process in each bin. These cross sections are, respectively, 
5.8 pb, 89.7 pb, 16.1 pb, 20.3 pb and 0.4 pb. 
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Figure 2: Distribution for data and Monte Carlo of a) the visible energy divided 
by the centre-of-mass energy; b) the sum of the transverse momenta of the two 
jets divided by the visible energy; c) the angle between the two jets in the plane 
transverse to the beams and d) the angle between the jets. The distributions in 
a) and b) refer to pre-selection level, those in c) and d) to the final selection. The 
arrows represent the position of the cuts, once all other pre-selection or selection 
cuts are applied. 
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Figure 3: Distribution of the visible mass for data and Monte Carlo after the appli- 
cation of all other selection cuts. 
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Figure 4: Measured cross sections for the various centre-of-mass energies, indicated 
by the points, compared with the Standard Model predictions, indicated by the 
band. The bars on the point show the sum in quadrature of the statistical and 
systematic uncertainties. The inner bars represent the statistical uncertainties. The 
width of the band corresponds to an uncertainty of 3% on the predictions, derived 
as discussed in the text. 
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